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ABSTRACT 



The object of this thesis is to develop a procedure 
for calculating the dynamic displacement response of 
the platform of a fixed offshore structure acted upon by 
a regular wave train. The structure considered has four 
legs in a square configuration, with the waves impinging 
normal to one side of the square. The procedure may be 
manipulated for use with other leg configurations and wave 
directions . 

The thesis is written in a manner useful to a designer 
of fixed offshore structures. An iterative procedure may 
be used to arrive at the critical wave -displacement 
combination. 

The types of waves considered are those for which Stokes ' 
3rd approximation applies. A modified version of Morison's 
theory is used for drag forces. Families of curves are pre- 
sented for use in predicting these forces and their centroids. 
The curves, plotted in dimensionless form, show the limits 
of applicability of Stokes* approximation. MacCamy- 

Fuchs diffraction theory is used for inertia forces. 

To calculate the dynamic response of the platform the 
authors fit the vibration problem of the structure to the 
classical theory for a linear, single -degree -of -freedom 
system. The wave forces are expressed in terms of a Fourier 
series. Displacements are calculated for twelve positions 
in one wave cycle, and a displacement curve is drawn. A 
simple tabular form is presented for calculating these dis- 
placements . 

Computer programs (In FORTRAN) are given for performing 
the bulk of the force and displacement calculations men- 
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tioned above. The theoretical procedure for unbraced and 
braced leg configurations is tested by experiment for 32 
cases. For these tests the ratio of maximum theoretical 
to experimental displacements varies from about 0.5 to 1.7. 

It is concluded that the designer should incorporate 
a factor of ignorance of 3 in bis design. 

The authors show that the largest waves are not 
necessarily the most critical for the structure. Rather, 
the critical design wave is a function of the frequency 
ratio (wave frequency/structure natural frequency). The 
designer must consider a range of significant waves which 
give frequency ratios of about 1.1 or less. 

It is recommended that calculations be performed for 
some real structures in ocean waves. From a comparison 
of the results with observed data a revised factor of ig- 
norance might be obtained. 

In Appendix I an approach to the solution for vortex 
shedding forces and frequencies in a train of surface 
waves is presented. Existing vortex theories are modified 
to apply to a vertical cylinder in finite amplitude waves, 
and boundaries are set on the predicted maximum forces . 
Experiments are conducted from which vortex forces are seen 
to fall within the boundary limits set by theory. The ex- 
periments also validate the proposed theory that the shedding 
frequency is a predictable integral multiple of wave fre- 
quency. It is recommended that a procedure be developed for 
determining platform displacements due to these forces. 
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I . INTRODUCTION 



The problem of wave -induced vibrations in fixed 
offshore structures suddenly was cast into the public 
news focus when on January 15, 1961, the United States 
Air Force Texas Tower Number Four collapsed 80 miles 
off the New Jersey coast in a heavy storm with the at- 
tendant loss of 28 lives [l3 • 

Prior to that incident, and since, no one has 
computed the dynamic displacements of such a structure 
acted upon by wave forces, to the best of the authors' 
knowledge. The design approach has been to make the 
natural frequency of the structure relatively high 
(about 50 cycles per minute ) and then simply treat the 
wave force as if it were statically applied. This 
method is crude at best, and can be very misleading. 

During the winter of 1958-59 a platform motion 
study [2] was made of Texas Tower No. 4. It was ob- 
served that the tower platform displacements were 
greater for 10 and ll ft. waves (about +2 in. ) than 
for 20 ft. waves. Obviously the highest wave was not 
the most critical for this structure. 

Numbers in brackets refer to literature citations. 
Appendix J. 
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The Air Force was interested in platform motions 
of Texas Tower No. 4 because they wanted to Irnow what 
errors might be introduced into 1 the installed aircraft - 
tracking radar system. (For this particular situation 
it was felt that rotational motion would probably be 
more of a problem than translational motion.) 

This thesis develops a method of calculating the 
dynamic translation displacements of a fixed offshore 
structure in a regular wave train. The theory developed 
has been experimentally checked by the authors. The 
experimental model was patterned somewhat similar to the 
United States Coast Guard's Buzzards Bay Light Station, 
located off the coast of Massachusetts £ 3 ] . 

The Coast Guard has initiated a long range pro- 
gram of replacing lightships by fixed offshore structures. 
The two completed to date (Buzzards Bay and Brenton Reef 
Light Stations) are braced structures with four legs. 

The U.S. Navy built a similar but considerably larger 
structure at Argus Island, near Bermuda. The U.S. Air 
Force built three Texas Towers, which were three-legged 
structures. Many other fixed offshore structures have 
been built in the past also, particularly in connection 
with drilling and mining operations in the Gulf of Mexico. 
The potential exists jj4^j for many more offshore oil- 
drilling rigs to be built around the world in the years 
to come. 
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There has been at least one case reported ^5^ 
in which vortex shedding was important in waves. In 
that case after a one week storm a two foot diameter 
pile suffered a fatigue failure caused by transverse 
vibrations of 2.5 seconds period in 12 feet high 
waves of 15 seconds period. The authors of this thesis 
suspected that vortex shedding forces might have a 
noticeable effect on offshore structures. Some investi- 
gation of this subject has been included in this thesis. 
The authors conceded from their literature survey that 
no previous work has been done on the question of vortex 
shedding forces on a stationary vertical cylinder in a 
wave train. 
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II . PROCEDURE 



A. Theoretical Procedure 
1 . Hydrodynamics 

General 

This section deals only with those forces on a 
cylindrical pile acting in a direction parallel to that 
of wave propagation. These forces are termed longitudinal 
forces to distinguish them from the transverse or "lift" 
forces caused by vortex shedding. The transverse forces, 
although requiring consideration in a practical structure 
design, caused such small displacements of the model under 
study that it was not feasible to obtain experimental 
verification of an analytic approach. Because of this, 
the transverse force investigation, both theoretical and 
experimental, is included as Appendix I and is not dis- 
cussed further in this section. 

The ultimate design procedure would be to make use 
of a statistical analysis of ocean wave and energy 
spectra in the development of the input forcing functions 
for the structure displacement problem. 
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Steps have been taken in this direction, for example 
[ 6 J, although the problems encountered, not only in ob- 
taining a sufficient scope of reliable data [ 7 J, but also 
in applying this data to the development of a useful de- 
sign criteria, have yet to be solved. For this reason 
the authors are restricting the investigation to a 
regular train of surface waves whose forces are capable 
of relatively simple mathematical superposition. 

In the past a typical design procedure has been to 
decide upon a "design wave," \fhich is the largest wave 
to be expected (within specified confidence limits) over 
the planned lifetime of the structure. Using applicable 
wave theory, one may then determine the wave forces and 
moments acting over an incremental length. These may be 
summed up in a tabular form to give the total force and 
moment. [ 8 , 3 ]. Some designs, also, have considered 
a breaking wave of size somewhat smaller than the design 
wave £ 9 , 10 J. 

If the dynamics of the problem are considered it 
is obvious that this design wave may not be the wave 
which causes the greatest excursions (or stresses) in 
the structure. Bather it is to be expected that a wave 
(of significant height) whose fundamental or some low 
order harmonic corresponds closely to the natural fre- 
quency of the structure may be the critical wave from 
a design standpoint. This possibility has been indicated 
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where, over an extended period of observation, a 7.5 foot 
wave caused greater platform motion of Texas Tower No. 2 



is cited in the Introduction for Texas Tower No. 4 . 

A serious situation also may exist when the wave 
length is approximately equal to the leg spacing, for 
in this case the fundamental wave forces (and possibly 
some harmonics) will be in phase. 

The proposed procedure, therefore, presents a 
method whereby several waves of various lengths and 
steepness may be considered over the range of possible 
water depths. The critical displacement -wave combination 
is then arrived at by an iterative process. To allow for 
a solution by this method within a reasonable length of 
time families of curves are developed such that total 
force and moment may be determined by a simple multi- 
plication. This is intended to supplant the tedious in- 
cremental tabulation and summation mentioned previously. 
The curves (Figures III, IV, V) are presented in terms 
of dimensionless parameters such that they may be applied 
to any situation for which the proposed theory is appli- 
cable . 

A pile spacing greater than 10 pile diameters has 
been found experimentally to be large enough such that 
proximity effects are negligible, [llj This spacing 
is exceeded, not only in the vast majority of offshore 



than much 




A similar observation 
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structures of the type considered, but also for the 
model which was the object of the authors' experi- 
mental study. 

Wave Theory 

H » wave height, measured from 
crest to trough. 

L * wave length. 

d = water depth. 

ryj Q m height of wave crest above 
still water level. 

Several theories on wave motions have been developed 
and tested experimentally with varying degrees of 
success. These theories are found to be applicable 
over finite ranges of wave parameters and water depths, 
with no one theory satisfying all conditions. Since 
the wav£s of interest in this study are of the re- 
latively shallow water finite amplitude type, two 
theories seem to be of particular interest - namely, 
Stokes' third approximation and the Solitary wave theory. 
These theories are discussed in detail elsewhere (for 
example £4, 12, 15 , 14]) and will not be repeated here. 
However, for purposes of evaluating this study, the 
following limits of applicability may be noted from 
these references. 
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Stokes' 3rd approximation: 



a. ) 


(7c/ H) 4 


0.625 


b.) 


(h/l) ^ 


0.142 


Solitary wave: 




a.) 


waves at or 


near breaking point 


b.) 


(d/L) — ^ 


0 


Theoretical limits on ocean wave parameters: 


a. ) 


(%/H) ^ 


0.75 


b.) 


(H/L) ^ 


0.142 


c. ) 


H„ = 0.78 
max 


(d-(H- ^ Q )) for shallow water 


d.) 


from the above, (H/d) ^ 0.653 



Using the maximum value of (/^? o /H) from Stokes' 
theory and the value of H^ above, one sees that 
(H/d) max is 0.603 by Stokes' theory. There is a 
gap between this value and O .653 for which the soli- 
tary wave theory would appear to be more applicable. 
However, waves have been observed breaking on Martha's 
Vineyard with values of H/d * O. 588 . [ 15 J 

Because this present investigation must be limited, 
it was decided to use only Stokes' 3rd approximation. For 

I 

the solitary wave and other theories the method of approach 
in determining forces and force centroids basically would 
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be the same, with the curve families of this thesis en- 
larged to include their particular ranges. Bretschneider 
£l4] has developed graphs which appear to be a good start- 
ing point for determining the range of applicability of 
the several theories, as well as values of H, L, and 
maximum crest velocity (U max ) for design waves. These 
values could serve as first trial inputs for the pro- 
cedure to be presented here. For a conservative estimate 
of the force exerted by a breaking wave on a piling a 
value of 4 times the force as determined from Stokes* 
theory for a wave of the same height may be used. [l6~j 

For reference purposes the applicable formulas 
from Stokes’ 3rd approximation are included here in 
dimensionless form. [ijJ The reader is referred to 
the list of symbols at the beginning of this thesis and 
Figure I for use with these equations. 

Let C » wave celerity, ft ./sec. 



J « 2rrd/L 

The vertical displacement of a particle from 
its mean position is given by (2) 



( 1 ) 




slnh 



(sinh J) 



i2J_S/dH 
nhJT J 



cos 29 



( 2 ) 
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Figure I . Definition Sketch of a Wave 



At the wave crest (2) becomes 



^max 

~a~ 



i $ sinh(J S/d) , 3 ir 
2^7 h[ ~£nh' J + H L/H 



sinh(2J S/d) 
(sinh J) 



and at the trough* 



(3) 



trough 

d 



JL_ f -sinh(J S/d) 3 jr sinh (2 J S/d) > 

-37H ^ sinh J + 5 i7H (slt * j W 



At the still water level (S/d - l), (2) gives the wave 
profile 
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3- rbf f 3ln 9 - [I !7h 002 zs j (5) 

From (5), with 9 ■ 90° the height of the crest above the 
still water level is 



%■- + Ii7h ^ 

and for 9 - 270° the distance to the trough is 





3 ir 

+ B I7H 



sinh(2J) } 
(sinh J) / 



The wave celerity is 



(6) 



(7) 



_C 

V gd 




J 7 1 + ( 



* \ 2 r 2 (cosh 2J) 2 + 2cosh 2J-tT7) /qa 

v** — si,)* — -if (8) 



The particle velocity at any depth S/d is 



u 

/gd 



X7H 




t 1 , 1 3 cosh (2 J S/d 

(sinh j) 2 2 (sinh J) 2 



)J 



cos 29 + 



jr 1 (* cosh (2J S/d) ^ 

iTif -2 / (sinh If JJ 



(9) 



At the wave crest (9) becomes 
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1 



u 

max 



IT C 




cosh(J S/d) 
sinh J + 



w 

175 * 



(slnh J)‘ 



1 .5 cosh(2J S/d) >> 
^ (sinh 7r ' 




cosh (2 J S/d) 
(sinh J r 




At the wave trough ( 9 ) becomes 



( 10 ) 



^trough ir C 
:32 — » T7w» 

jT#T V %/id 



cosh(J S/d) _tr 

sinh J I7 h* 



(sinh J)‘ 



1^3 cosh (2 J S/d) N 



(sinh J) 



+ 




( 11 ) 



Force Theory 

Using Stokes' 3rd approximation for wave motions, 
D.R.F. Harleman and W.C. Shapiro of the Massachusetts In- 
stitute of Technology Hydrodynamics Laboratory have de- 
veloped a force theory which shows good agreement with ex- 
periment. £ 17 J In developing this theory they applied 

certain modifications to Morison's Force Theory for drag 
and the MacCamy -Fuchs Diffraction Theory for inertia forces. 
In this thesis the authors will utilize the approach as pre- 
sented by Harleman and Shapiro. Families of curves will be 
developed from their equations, these curves to be used for 
the iterative design procedure mentioned previously. 
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As will be shown the only forces of interest are 
the maximum values of the drag and inertia components for 
both crest and trough regions. 

F D « maximum drag force in crest region 
(F D ) t * maximum drag force in trough region 
Fj « maximum inertia force in crest region 
(Fj ) t - maximum inertia force in trough region 
D » cylinder diameter 



Pi 



D 



Y D d 2 C D A 

* 2 



( 12 ) 



( F ) m -£2. d2 Cp trough B 

'u't 2 



( 15 ) 




1+V d 



where 



(«) 



o 



B 




o 



and S* is the vertical distance to the actual particle 
position and is therefore equal to the mean particle 
position S plus or minus the vertical displacement y. 



For the crest 



S' - S + y 



max 



and for the trough 



S' - S - y 



trough 



(16) 

(17) 



It can be seen that (14) and (15 ) do not lend themselves 
to regular Integration, but must be solved graphically or 
by some other means. The steady state drag coefficients 
C D and C D trough be used ln ( 12 ) and (13)* respectively, 
are determined as functions of Reynolds numbers, defined 
as follows: 



(R 



llF 

/ max ^ 


D J 


u D 

>> 


’ ^ V 



d A . 

V d 



n? 

v ; 



, R . v trough D , D , g . d B 

trough o ^ y 1 ~ ( %~ H ) 



(18) 

(19) 



P 



where u & and 

mcix 



u trough are r0 °^ mean squares 



of the velocity distributions for the crest and trough, 
respectively . 



J u 2 

V max 



A g d 
d + % 



( 20 ) 



J 



u 



trough 



B g d g 

d + (%-H) 



( 21 ) 



In (18) and (19) , A and B are determined from (14) 

t/j -H 

and (15)# respectively, and -y o /d, — from (6) and (7)* 
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Figure IX , Steady state drag coefficient versus 

Reynolds number for circular cylinders, 
fi7, iaj. 



respectively. Figure II is a plot of steady state vs, JR 
for circular cylinders. 



F 



I 



rf T D 2 H 
% cosh J 



sinh I cos 



Q 



( 22 ) 



2 

(F t ), * 77 ^ — T ' D t H sinh I cos 9 
' It 4 cosh J 

where I « J + sin 9 



(25) 

(24) 
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£17, 13J. 



respectively. Figure II is a plot of steady state vs. |R 
for circular cylinders. 

2 

F_ * ■ (, ^ y -r- D T - — sinii I cos © (22) 

I 4 cosh J . 

2 

(F t ). * i,^ H slnil 1 cos 6 (25) 

v I't 4 cosh J 

where I « J + -— 7 ^ sin © (24) 
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and 9 is that angle of wave position (0° < 9 <• 180° 



for (22), and 180 <■ 9 ^360 for ( 25 )) such that F-j. and 
(Fj )t are maximized, respectively. Because 9 will be 
ver*y nearly equal to 0° and l 80 ° for these situations, 
it can be seen from (24) that little error (and great 
simplification) is introduced by letting 9*0° and 
180 0 in (22) and (25) * respectively. 

This gives Fj - -(Fj^ - — ^ - T - Zf D - - H tanh J (25) 

Let S' «* center of action of the force F^ and 
(S’) t “ center of action of the force (F^)t 







1 + V d 



Then 



(26) 



o 



A 




o 



3 

Sj ■ center of action of force Fj 




sinh J - cosh 




(28) 
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Harleman and Shapiro have shown that S’ can be 
assumed constant for 0° ^ 9 ^ 180°, and (S' ) t can be 
assumed constant for 180° - 9 - 360° with little loss in 
accuracy. Also, with the assumption that S ^ is equal 
for crest and trough by (28), it is reasonable to assume 
that it is constant for any value of Q. 

By use of (5, 7, 10, 11,14,15,26,27) values of A and 
B for (12,13,18,19) and and for (26) and (27) 

were computed for families of the parameters H/L and d/L 
on an IBM 7090 digital computer. 

The results are plotted in Figures III, IV, and V 
in dimensionless form. Values of H/d equal to 0*603, 
which was mentioned as the approximate limit of applica- 
bility of Stokes' theory, and 0.653> which is the maximum 
value this parameter may obtain (page 8 ), are also 
shown in these figures. These values indicate the boundary 
for use of the figures, as well as the boundary up to which 
the solitary wave theory could be used to extend the range 
of applicability. The computer program, representative 
output, and comments are included as Appendix D. In ad- 

g | Q* | 

dltion to A, B, g- , and (g-) t > the following are given 
as outputs: 

% 

a - 

»-% 
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